Temporally ordered collective creep and dynamic transition in the 
charge-density-wave conductor NbSe 3 
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We have observed an unusual form of creep at low temperatures in the charge-density-wave (CDW) 
conductor NbSe3. This creep develops when CDW motion becomes limited by thermally-activated 
phase advance past individual impurities, demonstrating the importance of local pinning and related 
short-length-scale dynamics. Unlike in vortex lattices, elastic collective dynamics on longer length 
scales results in temporally ordered motion and a finite threshold field. A first-order dynamic phase 
transition from creep to high-velocity sliding produces "switching" in the velocity-field characteristic. 
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Interaction between internal degrees of freedom and 
disorder determines the dynamical properties of driven 
periodic media, a class of systems that includes vor- 
tex lattices in type II superconductors |Q, Wigner crys- 
tals g|, magnetic bubble arrays ||, and charge-density 
waves (CDWs), the low-temperature phase of quasi-one- 
dimensional conductors [Q . While the interplay between 
quenched disorder and elastic deformations is relatively 
well understood, a complete description of the role of 
thermal disorder and plastic deformations, which may re- 
sult in disordered dynamical phases such as driven glass, 
smectic and liquid states H , has not yet been achieved. 

CDWs have long been regarded as a prototypical sys- 
tem for the study of many-degree-of-freedom dynam- 
ics, both because of their relative theoretical simplicity 
and because CDW materials like NbSe3 exhibit collec- 
tive phenomena with remarkable clarity. A CDW con- 
sists of coupled modulations of the electronic density 
n = no + m cos[Q c a: + 4>{x)] and of the positions of the 
lattice ions [Q. Applied electric fields E greater than 
a threshold field Et cause the CDW to depin from im- 
purities and slide relative to the host lattice, resulting 
in a non-linear dc current density j c proportional to the 
CDWs sliding velocity. The impurities cause the CDW 
to move nonuniformly in both space and time, and the 
elastic collective dynamics leads to oscillations ( "narrow- 
band noise") in j c {t). The frequency v of these oscil- 
lations is proportional to the dc component of j c , and 
their Q — v/Av can exceed 30,000 in high-quality NbSe3 
crystals. 

Despite these simplifying features, most aspects of 
CDW transport at low temperatures remain poorly un- 
derstood. At temperatures T > 2Tp/3, j c above Et is 
a smooth, asymptotically linear function of E. However, 
at low temperatures j c {E) changes drastically, as illus- 
trated in Fig. |l|. CDW conduction still begins at Et but 
j c is small and freezes out with decreasing temperature 
for fields less than a second characteristic field E T > Et- 



At E T j c increases by several orders of magnitude to a 
more nearly temperature-independent value, often by an 
abrupt, hysteretic "switch." Similar behavior is observed 
in all widely studied CDW materials and is thus a fun- 
damental aspect of CDW dynamics fill 
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FIG. 1. (a) Form of j to t(E) in the CDW conductor 
NbSe3. Dotted line: single particle current density j 3 oc E. 
Dashed line: total current density jto-t, = 3" + 3e a ^ high 
temperatures (T > 2T P /3). Solid line: j tot at low T. The 
difference between the solid or dashed lines and the dotted 
line gives the CDW current density j c . (b) Temperature de- 
pendence of Et and Et in NbSe3. 

We have characterized the sliding CDW's transport 
and structural properties in the low-temperature regime 
of extremely high quality NbSe3 crystals. For Et < E < 
E T , we find that j c is activated in temperature and in- 
creases exponentially with field. Contrary to previous ob- 
servations in driven periodic media, this creep-like collec- 
tive motion exhibits temporal order. Our results illumi- 
nate the relation between local and collective pinning and 
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indicate that dynamics on lengths much shorter than the 
Fukuyama-Lee-Rice (FLR) length — neglected in most 
theoretical treatments — play a central role. They imply 
revised interpretations for "switching" at E T , the low- 
frequency dielectric response, low-field relaxation, and 
nearly every other aspect of the CDW response at low 
temperatures. 

High purity (t\r ~ 400) whisker-like NbSe3 single crys- 
tals with typical cross-sectional dimensions of ~3 by ~0.8 
/im were mounted on arrays of 2 /im wide gold-topped 
chromium wires []l0|| . The total current density j to t is a 
sum of the CDW and single-particle current densities, j c 
and j s . j c is orders of magnitude smaller than j s in the 
range Et < E < E T (except at relatively high temper- 
atures and very close to E T 0). j c (E) cannot be di- 
rectly measured and its form in the low-velocity branch 
of NbSe3 has not previously been determined. As shown 
in the inset to Figure 2, our high-quality crystals |ll[] ex- 
hibit voltage oscillations with Q's as large as 130 in this 
regime. Consequently, we are able to determine j c by 
measuring the oscillation frequency v = (Q c /27ren c )j c , 
where — e is the electronic charge and n c is the condensed 
carrier density. j c {E) was independently estimated by al- 
ternating the applied current's direction and measuring 
resistance transients R(t) associated with transients in 
the distribution of CDW strain e(x ) = (l/Q c )(dc/)/dx) 
between the current contacts |Io| , [l^] . 



fl o 

pi 

IP 



o 
O 



100 




a. 



o 



o 
o 



500 
field E 



700 
(mV/ cm) 



900 



FIG. 



sity 
Eq. 



i 



2. Coherent oscillation frequency v and current den- 
versus electric field E. The solid lines are a fit to 
The intersection of the lines with the horizontal axis 
corresponds roughly to the measured Et at each tempera- 
ture. The dotted vertical line indicates E T . Inset: spectral 
density S(f) at 22.8 K for E/E T = 2.63, 2.77 and 2.88; the 
curves are offset vertically for clarity. 



Figure shows the CDW current density j c {E) = 
v x 0.32 pA//im 2 Hz [[l3| calculated from the measured 
oscillation frequency v(E) for Et < E < E T at four tem- 



peratures. The CDW moves extremely slowly through- 
out this field and temperature range: the smallest mea- 
sured v values at T=20.7 K correspond to CDW motion 
of roughly one wavelength or 14 A per second and to 
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Between T«40 K and T«20 K, j c at fixed 



E < E T is temperature activated, decreasing by roughly 
7 orders of magnitude. j c jumps abruptly at E T , with 
j c (E=l.lE T )/j c {E = 0.9E T ) increasing from - 10 3 to ~ 
10 6 as T decreases from 28 to 20 K. 

The current density j c ex v can be fit by a modified 
form for thermal creep Jl4| 



j c {E,T) = a [E-E T ] exp 
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where the [E — Et] term describes the fact that the cur- 
rent drops to zero at a threshold Et that remains large 
even at high temperatures. The solid lines in Fig. ^| in- 
dicate a fit with To = 505 K, a = 136 KV _1 cm, and 
(To = 350 Sl^/im -1 . The value of To is insensitive to the 
assumed field dependence and corresponds to 0.6 times 
the single-particle gap 2A [O, consistent with measure- 
ments of delayed conduction |16|| and of a c near E T above 
30 K ||. Although creep is observed in other systems, 
the coherent oscillations imply that the creep in this case 
is highly unusual: it exhibits temporal order. 
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3. (a) Single-particle resistance R of a 70 /im seg- 



ment adjacent to a current contact versus electric field E. 
(b) R(t)/R(oo) for the same segment following a reversal 
of the polarity of E, as indicated by the arrow in (a), for 
E/Et=1A0, 1.69 and 1.93. (c) Comparison of j c calculated 
from R(t) |12j with j c obtained from measurements of the co- 
herent oscillation frequency. The current contacts were 630 
Hia apart, and £r(20.5 K) = 49 mV/cm. 

Figure || shows j c {E) at T=20.5 K obtained from 
transient measurements |U2|. These data agree closely 
with j c (E) deduced from v(E) over nearly three decades 
in j c p7[ |. Combining the two measurements yields 
jejv = 0.22 pA//im 2 Hz, consistent with the expected 
value of jc/v — 0.32 pA//zm 2 Hz |f| within the factor- 
of-two uncertainty in the value of j c determined from 
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transient measurements. This rules out significant fil- 
amentary conduction, observed in the low temperature 
regime of other CDW materials, and implies that the en- 
tire crystal cross section or at least a significant fraction 
of it is involved in coherent conduction. 

Figure [| shows the results of high-resolution x-ray 
diffraction measurements of the CDW's transverse struc- 
ture versus electric field. The CDW creates superlattice 
peaks in the diffraction pattern, and the half-width of 
each peak is inversely related to the CDW phase-phase 
correlation length. For E < Et, the resolution-corrected 
inverse half-width is I ~ 4100 A, comparable to the crys- 
tal dimension in this direction. For E > Et, I decreases 
monotonically with increasing E, remaining greater than 
2500 A for Et < E < E T . I does not show any abrupt 
change at E T despite the several orders-of-magnitude in- 
crease in j c there. Similar results were obtained in other 
directions perpendicular to Q c (e.g. [1 3] and [10 2]) 
and at higher temperatures. 
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FIG. 4. Inverse CDW peak half width (corrected for in- 
strumental resolution) in the [1 0] direction versus /tot- It 
and an upper bound for It were determined from measure- 
ments of dV/dI^ ^ and of the sharp increase in 1/f-like noise, 
respectively. 

Several different models have been proposed to ac- 
count for the low-temperature properties of CDW con- 
ductors. In K0.3M0O3 and TaS3, whose Fermi surfaces 
are completely gapped by CDW formation, the activa- 
tion energies for the single-particle conductivity a s and 
the CDW conductivity a c in the low-velocity branch are 
both comparable to the CDW gap so that <r c (T) oc cr s (T) 
M. Motivated by this observation, Littlewood Hj sug- 
gested that dissipation caused by single-particle screen- 
ing of CDW deformations limits CDW motion in the low- 
velocity branch, and that an abrupt, hysteretic transition 
to the high- velocity branch occurs at a frequency v com- 
parable to the dielectric relaxation frequency v\ oc a s 
when this screening becomes ineffective. The predicted 
value of v at the discontinuity for K0.3M0O3 and TaS3 is 
four orders of magnitude too large || , and for NbSe3 at 
T=20.7 K v is 13 orders of magnitude too large. Levy 
et al. |l(| showed that a related model exhibits a hys- 
teretic transition from the pinned state to a fast sliding 
state when a s is small even if high-frequency screening 
effects are neglected. Neither model can explain the low- 
temperature CDW properties of partially-gapped NbSe 3 , 



for which a s remains metallic and increases with decreas- 
ing temperature below 50 K. 

Various forms of CDW plasticity including phase 
slip at isolated defects flpfl and shear between two- 
dimensional CDW sheets 119fl have been suggested to ac- 
count for the properties of the low-velocity branch and 
the transition at E T . Our observation of highly coherent 
oscillations in high-quality crystals and earlier results j|] 
rule out models based on slip at rare isolated defects and 
contacts, and our x-ray measurements rule out the form 
of shear plasticity discussed in Ref. Q . 

Brazovskii and Larkin |^0| have focused on the CDW's 
local interaction with defects. At low temperatures CDW 
phase advance past rare defects occurs via thermally- 
activated soliton generation, and motion becomes much 
more rapid at large fields when the effective barrier to 
soliton generation vanishes. This interpretation has ap- 
pealing features, but the suggested form for the j c (E) 
relation at low temperatures does not reproduce the two 
branches separated by an abrupt hysteretic transition or 
the field dependence in either branch observed experi- 
mentally in NbSe3. Furthermore, the predicted E T is 
determined by the soliton energy and should be indepen- 
dent of crystal size. Experimentally, in NbSe3 both Et 
and E T vary as 1/t for crystal thicknesses t less than ~ 
20 /im ||. The thickness dependence of Et results be- 
cause transverse CDW correlations are limited by t so 
that collective pinning is two-dimensional [ pT| . Conse- 
quently, the thickness dependence of E T implies that it, 
too, is determined by collective effects. 

CDW creep of a fundamentally different character is 
observed in thin NbSe3 crystals at high temperatures 
pl| , p2| . Near Tp, Et is rounded, nonlinear conduction 
can extend to near E=0, and highly coherent oscillations 
below the nominal Et are not observed. This incoherent 
creep occurs when kpT approaches the collective pin- 
ning energy (oc A(T)H) of the phase-correlated FLR 
domains, which in NbSe3 have micrometer dimensions. 
The temporally-ordered creep observed in relatively thick 
crystals at low temperatures above a sharp threshold Et 
must involve barriers that are much smaller than those 
of collective pinning and that are not rare, and a length 
scale that is much smaller than that of the collective dy- 
namics responsible for the narrow-band noise. 

Motivated by earlier ideas J20]j2^ , |24]| , we suggest that 
the low-velocity branch develops when CDW motion be- 
comes limited by thermally-activated phase advance by 
~ 27T past individual impurities. Although collective pin- 
ning is weak [ElJ, the phase of the Q c = 2kp oscillations 
is fixed at each impurity so that phase advance requires 
CDW amplitude collapse and a finite barrier ~ A [ pi[ . 
Collective dynamics within volumes containing enormous 
numbers of impurities (set by the FLR length) then gen- 
erates the finite threshold Et and coherent oscillations, 
as in the high-temperature regime. Unlike in vortex lat- 
tices, long-length-scale CDW dynamics is largely elastic 
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and thus retains temporal order even though the short- 
length-scale dynamics is stochastic. 

The E — Et prefactor J2j| and the remaining terms 
in Eq. |l| follow naturally from this combination of long 
and short length-scale processes. The measured barrier 
To is consistent with the expected pinning barrier per 
impurity of ~ A p3]. An applied electric field should 
reduce this barrier by ~ en c VXE, and using a conden- 
sate density n c = 2x 10 21 citT 3 ppl|, a CDW wave- 
length A = 14 A and the measured a value yields a 
volume V involved in each thermally-activated event of 

V f» 4.2 x 10~ 17 cm 3 . Using the scale factor expected 
for typical impurities (2^j2^] , the bulk residual resistance 
ratio of ~ 400 for our crystals corresponds to a concen- 
tration rii w 2.5 x 10 16 cm~ 3 . The volume per impurity 
1/rii w4x 10~ 17 cm 3 is thus in excellent agreement with 

V deduced from creep measurements. 

The present experiments together with those of Rcf. ||] 
rule out all previous explanations of the "switching" be- 
tween low and high velocity branches at E^, in NbSe3. 
We suggest that switching occurs via a first-order dy- 
namic phase transition ||^7||. The long-length-scale dy- 
namics exhibits temporal order in both branches, but in 
the high-velocity branch dynamic fluctuations produced 
as the CDW moves past impurities may become more 
important than thermal fluctuations in overcoming impu- 
rity barriers p7j |. The transition's abruptness, hysteresis, 
and temperature dependence shown in Fig. [j], the CDWs 
tendency to fragment near the transition into distinct 
conducting regions and the field and temperature- 
dependent time delays required for the transition's com- 
pletion Jl6| are all consistent with this explanation. 

Finally, we note that local temporal order has very 
recently been observed in the creep regime of a vortex 
lattice |||. 
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